Relativistic Ferromagnetic Magnon at the Zigzag Edge of Graphene 
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We study the spin-wave excitations near the zigzag edge of graphene. It is rather interesting 
that we obtain a single branch of relativistic ferromagnetic magnon due to the presence of the open 
boundary. Note that magnons in antiferomagnets appear in pairs, while the single brach magnon 
in ferromagnets does not have relativistic dispersion. Thus, the magnon near the zigzag edge of 
graphene is a hybrid of both, signaling its intrinsic property as a boundary excitation that must be 
embedded in a higher dimensional bulk system. 
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Graphene is the two-dimensional single-layer graphite 
composed of carbon atoms in honeycomb lattice. Re- 
cently, this novel material was successfully fabricated in 
laboratories [U and stimulates intense investigations 
on both experimental and theoretical sides. Its surpris- 
ing band structure and high mobility seem promising for 
potential applications in many aspects[4J. For instance, 
there are proposals 5 6 that graphene nanoribbons can 
be the building block for future quantum computation. 
Since graphene is a low dimensional system, we expect 
the electron-electron interaction should play some role. 
In a recent paper, Son, Cohen and Louie [7] spotted local 
magnetic moments near the zigzag edges of a graphene 
nanoribbon due to the Coulomb interaction. In fact, 
correlation-induced magnetic moment in graphene has 
be speculated for quite a while [5J H EH ED III] but its 
low-energy excitations remain poorly understood at this 
point. 

The novel magnetism shows that electronic correla- 
tions play a significant role at the edges of graphene. In 
condensed matter systems, the physical properties at the 
edge are often tied up with related bulk properties. For 
instance, the Andreev bound state E2 near the edge of a 
superconductor are related to the edge topology and also 
the pairing symmetry in the bulk. The edge states pH] of 
a quantum Hall liquid are described by the chiral Lut- 
tinger liquid that cannot be separated from the bulk due 
to gauge invariance. Therefore, we are inspired to inves- 
tigate the low-energy spin excitations near the edge of a 
single-layer graphene here. For simplicity, we concentrate 
on a semi-infinite honeycomb lattice with zigzag edge as 
shown in Fig. [I] Since the band structure of graphene is 
well approximated by the nearest-neighbor hopping t, the 
tight- binding model is sufficient. Furthermore, if we as- 
sume the electron-electron interaction can be capture by 
an effective on-site interaction U, it is natural to model 
the semi-infinite graphene by the Hubbard model at half 
filling (n) = 1. 

Unfortunately, it is known that the two dimensional 
Hubbard model cannot be solved exactly/reliably at cur- 
rent stage. Thus, approximations are inevitable. Since 
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FIG. 1: Lattice structure for the semi-inifinite graphene with 
zigzag edge. The Neel state is shown by the alternative spins 
on different sublattices. 



we are interested in the spin physics, we propose to con- 
centrate on the effective Heisenberg exchange coupling in 
the system, 



H= J^S(1).S(1'), 

<i,i'> 



(1) 



where 1 labels all lattice sites and (1, 1') means nearest 
neighbor pairs. The coupling constant is J = 4t 2 /U > 
in the strong interaction limit but is viewed as an in- 
dependent parameter here. The honeycomb lattice of 
graphene can be separated into two sublattices A and B. 
The outmost edge sites belongs to the sublattice A in 
our convention as shown in Fig. [T| Starting from a Neel 
ordering state, it is convenient to rotate the spin on sub- 



lattice B by S x — ► S x , S y 



-S y , and S z 



-S, so that 



the semi-classical ground state become uniformly ferro- 
magnetic. Following the standard procedures, we can 
represent the spin on the lattice sites by the Holstein- 
Primakov boson 



S+(r) = y/2S-b[(r)b L (r) b L (r), 

SZ(t) = b{(r) y/2S-bl(r)b L (T), 

S£(r) = S-b[(r)b L (r), (2) 

where &z,(r), b L (r) are annihilation/creation operators for 



bosons. To make the sublattice dependence explicit, we 
have separate the lattice positions into sublattice index 
L = A,B and r for the triangular lattice. Expanding 
the Heisenberg interaction to the quadratic order in the 
bosonic operator, one can obtain the familiar Hamilto- 
nian for spin- wave excitations, 

H = JSY, b A (r)b A (r) + b B (r + X)b B (r + A) 

r,A L 



+b A (v)b B (r + X) + b B (v + X)b A (v) 



(3) 



The summation over A includes all three nearest neigh- 
bors \i = a(±l/2, l/2\/3),a(0, — 1/V3) as shown in 
Fig. [l] The presence of open boundary ruins the trans- 
lational invariance in y direction and the above Hamilto- 
nian cannot be diagonalized by the usual Fourier trans- 
formation. 

However, since the translational invariance in the x 
direction is still valid, we can perform the partial Fourier 
transformation 



1 



(4) 



to simplify the Hamiltonian. After some algebra, the 
spin-wave Hamiltonian can be casted into the form, 



H 



h-21 h-22 



(5) 



where we introduce the two-component spinor ^ = 
[if> A (k x , y), ip B (— k x , y)]- The matrix elements in the re- 
duced 2x2 matrix are semi-infinite matrices: hn — 



2 + £>t n h 



22 



3 and h 2 i 



J i + J 2D, where 



Ji = 2 cos %^ comes from the tilted bonds and J2 = 1 
from the vertical bonds. For anisotropic interaction, J2 
would derivate from one. The semi-infinite displacement 
matrix D is 



10 

10 

D= 1 





Now the problem in converted into solving the eigen- 
states for a one dimensional semi-infinite system. For no- 
tation clarity, we suppress the k x momentum dependence 
in the following. Writing down the matrix elements of 
the Hamiltonian explicitly, it leads to the coupled Harper 
equations, 

3^u(ra) + Jiip B (n) + J2ipB(n - 1) = eip A (n) 
-3ip B (ri) - JiipA(ri) - J2ipA(n+ 1) = eip B (n), (6) 

where n = 1, 2, 3, ... denote the lattice coordinates in the 
y direction. However, one needs to pay special attentions 




FIG. 2: Dispersion for bulk and edge magnons with JS = 1 
and a — 1 simplicity. The shaded area represents the bulk 
magnons with double degeneracy AS Z = ±1. The thick blue 
line indicates the dispersion for the edge magnon with only 
one brach AS Z = — 1. 



to the boundary. For zigzag edge, there is one missing 
bond from the outmost A sites to the opposite sublattice. 
The semi-infnite nature of the Hamiltonian gives rise to 
the constraint, 



iki(l) + JalMQ) = 0. 



(7) 



Now we are ready to find out the eigenstates by the 
generalized Bloch theorem |20j. Note that eigenvalue of a 
unitary displacement operator can be written as z = e , 
where k is the momentum. However, the presence of open 
boundary changes to this nice connection to the momen- 
tum. It is easy to check that the displacement operator 
satisfies DD^ = 1 but not in the inverse order D^D^l. 
Thus, it is no longer unitary. Thus, the eigenstate should 
be written in the more general form, 



* = 







CA 






c B 



(8) 



Upon substitution into the coupled Harper equations, the 
eigenvalue z of the displacement operator satisfy the sim- 
ple algebraic constraint, 



J1J2 ( z + 



1 



+ (J? + Jf + e 2 - 9) = 0. 



(9) 



Usually, for given momentum k, we look for energy e 
of the system. Here it is more convenient to look for 
the complex z for a given energy e. It is obvious that 
the solutions z, l/z appear in pairs. Thus, the solutions 
can be classified into two types: (1) plane-wave solutions 
with z — e ±zk (2) evanescent modes with real z — a, 1/a. 
Since we require the wave function to be finite at infinity, 
only one of the evanescent solutions survives while, not 
surprisingly, both plane- wave solutions are allowed. 

The plane-wave solutions correspond to the bulk 
magnons. The presence of the open boundary does cause 
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FIG. 3: Spin density profile for (a) k x — 0.1 and (b) k x — 0.5 
edge magnons. For clarity, the magnitudes on the edge are 
held constant. The magnitudes of the spin density (after sub- 
tracting the background from the ground state) are repre- 
sented by the radius and also the color gradients. 



a rather complicated mixing between counter propagat- 
ing modes ±k y . However, the energy spectrum remains 
the same as shown in Fig. [2] Therefore, we would only 
concentrate on the more interesting evanescent modes 
where z is real. Solving the coupled Harper equations 
together with the boundary condition, there exist only 
one solution e = \J 4 — Jf. Therefore, the dispersion for 
the single-branch edge magnon is 



E e (k x ) = 2JS\sm(k x a/2)\, 



(10) 



which becomes linear (rclativistic) in the small momen- 
tum regime. Note that, as shown in Fig. [2j the velocity 
of the edge magnon is slightly smaller than that for the 
bulk magnon and thus is protected by a gap. In fact, 
this protection arises from the sharp distinction between 
\z\ = 1 and \z\ < 1 solutions. Furthermore, the wave 
function of the edge magnon is 



1/z 
-1 



with 



- | sm(k x a)\ 
cos(k x a/2) 



(11) 



Note that < z < 1 as required for the evanescent modes 
except at k x = where the edge magnon merges into the 
bulk with z = 1. 

The spin density profile of the edge magnon can be 
obtained by subtracting the background form the ground 
state, AS* 2 (r) = (l|S z (r)|l) - (G\S z (v)\G), where |1) is 
one magnon state and \G) is the ground state. Since the 
spin density is uniform in the x direction, a factor l/N x 
appears where N x is the total transverse length. The 
spin density can be expressed in terms of the normalized 
wave functions, 



AS z (r) 



"ifl^AHl 2 , re A, 
+ ^\ip B (n)\ 2 , veB. 



(12) 



Summing over all lattice sites, it is easy to check the 
edge magnon we found carries the total spin AS Z = 



^ r AS z (r) = —1, as for those found in conventional fer- 
romagnets. The spin density profiles of the edge magnon 
with two different momenta k x — 0.1 and k x = 0.5 are 
shown in Fig. [3j For momentum fairly close to k x = 0, 
the spin density extends into the bulk with long decay 
length. However, for k x — 0.5, though not very far away 
from the origin, the spin profile is already quite localized 
to the zigzag edge. 

So, what's the significance of the above simple calcula- 
tions? For an ordinary ferromagnet, there is one branch 
of non-relativistic magnon. On the other hand, antiferro- 
magnetic ground state supports two branches of relativis- 
tic magnons. Here, near the zigzag edge of graphene, we 
find a hybrid between the two cases - a single branch of 
the relativistic ferromagnetic magnon. This implies that 
the spin-wave excitations near the edge of the graphene 
cannot be described as the ordinary one dimensional fer- 
romagnet and shall be viewed as an intrinsic boundary 
field theory. The missing branch can be understood by 
the following symmetry argument. Although an antifcr- 
romagnetic order breaks both time-reversal and sublat- 
tice exchange symmetries separately, the ground state 
remains invariant under both transformations simultane- 
ously. Thus, starting from the AS Z = —1 branch, one 
can reverse the spin orientations by time-reversal trans- 
formation followed by sublattice exchange to construct 
the other AS Z — 1 branch. This explains why antiferro- 
magnetic magnons on bipartite lattices always appear in 
pairs with opposite spins. However, the sublattice sym- 
metry is broken in the presence of the zigzag edge, this 
makes the single-branch magnon possible. In fact, we 
check with the armchair edge which doesn't destroy the 
sublattice symmetry and found no evidence for the novel 
edge magnon discussed here. 

Finally, we address the validity of the approximations 
in the linearized spin-wave theory. First of all, the charge 
excitations are all left out. Since there is no charge gap 
in graphene, it is expected that the neutral particle-hole 
excitations with AS Z = —1 will mix up with the edge 
magnon. However, since the Fermi surface shrinks to 
two Dirac points, the phase space is largely reduced. Fur- 
thermore, since the Fermi velocity vf is expected to be 
faster than the magnon velocity v m , the hybridization is 
further suppressed by a finite gap A ~ (vp — v m )k x for 
any finite momentum. This is similar to the protection 
between bulk and edge magnons we discussed previously. 
Another approximation is the existence of the Neel or- 
der. This is the most serious approximation since it is 
likely that the antiferromagnetic order may not exist at 
all in graphene. Though the quadratic fluctuations from 
the magnons do not destroy the Neel order and the spin- 
wave theory remain self-consistent, it remains an open 
question whether our findings here are true beyond spin- 
wave theory. Given the novelty of the edge magnetism 
in graphene, alternative theoretical investigations are in 
oder. And, as always, more experimental data will help 



to clear up our understanding of spin-wave excitations 
living on the edge. 
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